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INTRODUCTION 


Waves  on  the  ocean  surface  exemplify  a  number  of  important  lessons 

in  the  mathematical  study  of  random  processes.  The  casual  observer 

sees  a  series  of  irregular  moving  crests  and  troughs  that  gradually 

grow  and  shrink  in  the  course  of  time.  Mathematically,  the  ocean 

surface  within  certain  boundaries  is  a  single-valued  function  of  the 

planar  coordinates  (y,z)  and  time  (t) .  Let  it  be  X(y,z,t).  The  air 

motions  above  this  surface  and  the  water  motions  below  it  must  satisfy 

the  nonlinear  boundary  conditions  at  the  surface.  These  equations 

have  never  been  solved  exactly.1  The  reasonable  approximations  that 

2 

yield  solutions  corresponding  to  experiences  have  been  determined. 

Yet  even  under  the  simplifiying  assumptions  that  yield  provisional 
solutions,  the  behavior  of  X(y,z,t)  indicates  a  complex  structure 
involving  distances  ranging  from  centimeters  to  many  kilometers. 

A  one-dimensional  wave  study  begins  with  a  continuous  recording  of 
the  sea  surface  elevation  at  a  particular  point  (y 0,zQ) ,  the  recording 
instrument  being  functionally  equivalent  to  a  "blind  staff"  which 
ignores  the  directions  from  which  the  waves  approach.  If  an  electro¬ 
mechanical  pressure  transducer  were  placed  on  the  bottom  at 
(x0,y0,zQ),  one  can  define  xo»0,  and  the  instantaneous  pressure  would 
be  linearly  related  to  the  varying  depth,  which  is  X(y0,zQ,t).  Let 
V(t)  be  the  voltage  out  of  the  transducer.  Then 


t 

V(t)  =  /  X(y0,zo,t')h(t-t'jM)dt' 


(1) 
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expresses  V(t)  as  the  convolution  of  X  with  h,  the  impulse  response  of 
a  low-pass  filter  that  models  the  interaction  of  the  surface  with  the 
device  at  average  depth  \i.  Very  long  term  observation  of  V(t)  shows 
that  it  has  an  average  value,  proportional  to  the  mean  depth,  to  which 
periodic  and  quasi-periodic  fluctuations  are  added.  In  practice,  the 
duration  (Tr)  of  the  recorded  waveform  will  be  insufficient  to 
define  the  true  mean.  For  example,  a  15  minute  record  will  yield  a 
value 


v  =  (1/Tr)  f  V(t)dt  (2 
0 

which  reflects  the  short-term  average  depth  for  a  given  phase  in  the 
tidal  cycle.  It  makes  sense  then  for  the  analyst,  who  is  interested 
in  the  activity  of  the  waves,  to  take 


X  (t)  *  [V(t)  -  v]/c 

as  the  history  of  sea  surface  fluctuation  at  the  point  in  question, 
where  c  is  a  constant  that  accounts  for  the  midband  behavior  of  the 
low-pass  filter  model.  In  other  words,  X(t)  is  the  one-dimensional 
fluctuation  as  a  function  of  time. 


OCEAN  WAVE  STATISTICS 


The  two  most  salient  features  of  any  wave-like  phenomenon  are  its 
amplitude  and  period.  For  a  pure  (sinusoidal)  wave, 

Asin(2nt/T  +  0) , 

the  amplitude  (A)  ,  period  (T) ,  and  phase  (P)  can  be  determined  by 
inspection  of  a  single  wavelet,  as  the  behavior  of  X  on  (tjyt^+T) 
is  repeated  during  (t^+T,tj+2T) ,  and  so  on.  Ocean  waves,  however. 


WWW 
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are  only  quasi -per iodic;  and  the  peak  fluctuation  at  successive  crests 
is  likewise  subject  to  variation. 

Based  on  physical  insight  into  the  generation  and  propagation  of 
ocean  waves,  or  on  an  understanding  of  abstract  harmonic  analysis,  the 
scientific  observer  knows  that  X(t)  can  be  represented  very  precisely 
over  any  observation  interval  (of  Tr  seconds)  as  a  sum  of  sinusoidal 
components: 


X(t)  =  £  Ansin(2irnt/Tr  +  0n)  (4) 

n=l 

Since  the  flow  of  energy  that  the  component  wave  embodies  is 
proportional  to  the  squared  amplitude,  the  sequence  {An,  n-1,2,...} 
shows  the  allocation  of  total  power  among  the  components.  The 
mathematical  tools  which  extract  the  component  amplitudes  from  the 
recordings  of  X(t)  fall  into  the  domain  of  spectral  analysis.  Fast 
Fourier  Transform  (FFT)  algorithms  are  commonly  used  today. 

Most  of  the  significant  studies  of  ocean  wave  statistics  were 
written  prior  to  the  widespread  availability  of  the  digital  electronic 
technology  that  makes  FFT  the  standard  analytic  tool  that  it  has 
become.  These  pre-1965  studies  often  assume  a  tautological  model  in 
which  X (t)  is  represented  as  a  sequence  of  individual  wavelets.  Let 
t^  be  the  time  of  the  k-th  up-crossing  of  the  axis  (X*0)  by  X(t). 

Then  if  the  recording  begins  with  an  up-crossing  at  t=0,  the  sequence 

fc2 '  t3_t2'  t4~t3 '  **•'  tK"tK-l 

gives  the  respective  periods  of  the  successive  wavelets.  Defining 
T^  *  the  tautological  model  is 

K 

X(t)  -  23  B|<sin  (2*t/Tk) 1  (fck-l— fc<tk^ 
k-1 


(5) 
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where  the  indicator  function  l(*)  is  one  when  its  argument  condition 
is  satisfied,  and  zero  otherwise.  To  determine  B^,  one  could  take 
^ 2  times  the  root-mean-square  (r.m.s.)  fluctuation  over  the  corre¬ 
sponding  interval,  or  just  the  peak  positive  fluctuations,  or  some 
other  measure.  Indeed,  there  is  no  way  to  choose  the  B^'s  that  will 
make  equation  (5)  exact,  since  the  sine  function  is  just  an  approxi¬ 
mation  of  the  behavior  of  the  ocean  wave  over  each  period.  While 
equation  (4),  on  the  other  hand,  may  be  precise  in  principle,  the 
number  of  terms  in  the  summation  will  in  practice  be  truncated  at  some 
value  N: 


AjjSin  (2irnt/Tr  +  0n)  (4* 

example,  include  only  the  N  most  significant 

N 

*  Sn 

n*l 

adding  up  to  90  percent  of  the  total  energy  in  the  record.  Thus, 
whether  the  analyst  uses  the  tautological  model  of  equation  (5)  or  the 
spectral  linear  model  of  equation  (4'),  the  model  is  just  an  approxi¬ 
mation  to  the  true  behavior  of  X(t). 

Bretschneider  observes  that,  for  application  to  the  design  of 
coastal  engineering  structures,  there  is  a  greater  need  for  knowledge 
of  the  joint  probability  distribution  of  wave  heights  and  periods  than 
for  the  actual  wave  spectrum.**  in  this  and  similar  applications, 
the  effects  of  waves  are  nonlinear  functions  of  size  and  duration. 
Hence  the  tautological  model  is  preferred.  It  is  usually  assumed  (and 
generally  borne  out  by  analysis)  that  B^  and  T^  are  independent 
random  variables.  In  other  words,  the  (marginal)  distribution  of  the 
random  variable  Bk  is  unaffected  by  knowledge  of  Tk,  and  conversely 


N 

x(t)  i  y 


The  analyst  may,  for 
components,  with 
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Since  (T^ ,  k=l,...,K>  is  a  sequence  of  independent,  identically 
distributed  random  variables,  it  is  important  to  know  the  probability 
density  function  (p.d.f.)  of  T,  a  randomly  selected  element  of  the 
sequence.  To  this  end  one  defines 


t  =  T/T 


where 


T  =  ZTk/K  <7 

k=l 

is  the  mean  period.  The  variance  of  t  will  then  reflect  the 
deviation  from  pure  sinusoidal  behavior.  Furthermore,  it  has  long 
been  known  that  the  shape  of  the  wave  period  distribution  is  somewhat 
more  consistent  from  place-to-place  and  time-to-time  than  is  its 
centroid.  Therefore,  the  distribution  or  p.d.f.  of  the  relative  wave 
period,  t,  has  broader  generality.  Many  experiments  through  the 
years  have  shown  that  the  p.d.f.  of  the  relative  wave  period  is 
approximated  by 


p(x)  =  KjLT3e"blT  (8 

2  2  2 

where  K-^  *  ir/Ex  ,  bj^  =  tt/4Et  ,  and  Ex  «  1.07871.  The  first,  second, 
and  third  quartile  points  of  this  distribution  are  0.8,  1.0,  and  1.2, 
respectively,  so  that  half  of  all  the  periods  are  within  +20%  of  the 
median. 

The  definition  f  =  1/T  of  the  "frequency"  of  a  wavelet  leads  to 
consideration  of  the  relative  frequency  v  *  f/f  ■  1/x,  by  analogy 
to  equation  (6).  It  follows  from  equation  (8)  that  the  p.d.f.  of  v 


Pv<s) 


K2S-5e-b2/s« 


•  -  -j-  v.vV  / 


s*.  /aW  nWa’ 
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This  is  the  Bretschneider  spectrum  of  ocean  waves.  When  (9)  applies, 
half  the  energy  of  the  waves  is  concentrated  at  frequencies  between 
0.75f  and  l.lf.  Yet  the  relationship  of  the  Bretschneider  spectrum  to 
the  Fourier  spectrum  of  ocean  waves  is  not  a  simple  matter.  The 
transformation  which  converts  (8)  to  (9)  is  consistent  with  the 
tautological  model,  but  not  with  the  linear  model.  As  an  example, 
consider  the  case  of  white  Gaussian  noise  for  which  the  relative 
period  has  the  p.d.f. 

,  V  -T 

w  ( x )  =  e 

_2 

The  change  of  variable  v  =  1/t,  Idtl  =  v  dv,  gives  the  corresponding 
v- spectrum 

wv(s)  =  s“2e"1/s, 

which  implies  that  half  the  energy  is  distributed  between  relative 
frequencies  of  0.7  and  3.5.  But  the  Fourier  spectrum  of  white 
Gaussian  noise  is  flat,  since  a  sufficiently  long  record  of  this 
process  is  delta-correlated  and  the  Fourier  transform  of  the  delta 
function  is  a  constant.  On  the  other  hand,  inverse  Fourier 
transformation  of  the  last  expression,  taken  as  a  spectral  density 
with  s  =  2nif/f,  would  give  an  autocorrelation  function  of  the  form 
(1/fu)  Kj.  (2u)  ,  where  u  * ‘^iriFt  and  K.^  is  the  modified  Hankel 
function  of  first  order. 

The  Bretschneider  spectrum  is  only  one  of  the  theoretical  ocean 
wave  spectra  that  have  been  formulated  over  the  years.  Walden 
portrays  a  unanimous  disagreement  among  the  different  theoretical  wave 
spectra,  remarking  that  "it  is  strange  that  the  authors  of  the  various 
hindcasting  methods  nearly  always  find  that  their  data  agree  well  with 
the  properties  of  the  sea  in  question  even  though  the  differences 

_4 

between  the  diagrams  or  spectra  are  immense. 
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Then,  on  the  X  -  Y  plane,  the  joint  density  function  can  be  written  in 
terms  of  polar  coordinates  r  and  6: 


2  2 

Px (x)Py (y)dxdy  =  (2ira2)-ie“r  /2o  rdrd0  . 

Integrating  through  6  on  the  right,  one  obtains 

2  2 

pr(r)  =  (r/o2)e“r  /2o  ,  (10 

the  Rayleigh  density. 

The  same  result  can  be  derived  using  the  more  familiar  techniques 
of  probability  theory  and  real  analysis.  With  reference  to  the 
tautological  model  of  equation  (5) ,  one  has 

X  =  B  sin  (a>t) 

for  every  wavelet,  and 

Pr{  IXI  <xlB*b}  =  (2/ir) arcsin  (x/b)  (11 

when  t  is  uniform  on  the  interval.  Since  I X l  is  the  symmetric  and 
B  is  non-negative  by  convention. 


Pr{ IX  I <x}  *  f  Pr{ lXl£xlB*b}f (b)db, 
o 

where  f (b)  is  the  (unknown)  p.d.f.  of  the  amplitude  (B) .  Breaking  the 
integral  in  two  about  the  point  x,  the  last  equation  is  the  same  as 

x  00 

Pr{lXl<x>  *  /  Pr{ IX I <xlB=b)f  (b)db  +  /  Pr{ I X I <xlB«b}f (b)db. 
o  x 

But  IX l  is  clearly  less  than  x  of  b  is  less  than  x.  Therefore, 
the  probability  measure  in  the  first  integral  is  unity.  Hence, 


*  ■>  V-  _*  _»  -> 


CM 
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Pr{lXI<x)  *  £  f (b)db  +  £  2  arcsin (j£) f (b) db 


where  equation  (11)  has  been  substituted  in  the  second  integral.  The 
p.d.f.  of  IX I  is 

g  (x)  =  <L_  Pr  { I X  l  <x} 
dx 

by  definition.  It  follows  from  differentiation  of  equation  (12)  and 
an  integration-by-parts  that 


g  (x)  =  f  — 2f  (b)  db 

x  n  V  b2  -  x2 

The  question  now  becomes,  "If 


V”  2  2 

I  O  e~X  /2a  •  x  -  °' 


then  what  p.d.f.  f(b)  jointly  satisfies  (13)  and  (14)?"  The  answer  is 
given  in  equation  (10),  i.e.. 


f(b)  =  pr(b)  =  (b/o2)e-t>2/2o2  . 

The  derivation  follows  from  the  change-of-var iable 


a2  *  b2  -  x2 ,  db  =  ada/  Va2  +  x2  ; 


but  the  details  are  omitted  for  the  sake  of  brevity. 

The  equivalence  of  narrowband  Gaussian  fluctuations  (on  the  one 
hand)  to  a  Rayleigh-modulated  sinusoid  (on  the  other)  is  consistent 
with  both  the  linear  model  and  the  tautological  model.  Yet  it  should 
be  kept  in  mind  that  this  equivalence  is  contingent  in  either  case 
upon  the  spectral  width  being  small  as  a  fraction  of  center 


(10’) 


JN 


■>  ■>  V  V.’-‘  *-• 
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frequency.  This  proviso  is  not  always  satisfied  by  ocean  waves. 
Cartwright  and  Longuet-Higgins  have  derived  a  family  of  p.d.f.'s  to 
describe  the  distribution  of  ocean  wave  amplitudes.  The  Cartwright- 
Longuet-Higgins  density  is  parameterized  by  an  e  which  measures  the 
relative  spectral  width  of  the  fluctuations.  In  the  case  of  e*0, 
corresponding  to  a  very  narrow  spectrum,  the  C-L-H  density  coincides 
with  the  Rayleigh.  As  e  approaches  one,  the  p.d.f.  assumes  the 
Gaussian  shape  with  zero  mean.  In  this  latter  case,  the  envelope  of 
the  fluctuations  crosses  zero  from  above  and  below  in  succession,  as 
illustrated  by  Figure  2. 


DATA  ANALYSIS 


Figures  3  through  22  show  some  ocean  wave  records  together  with 
spectra  and  normalized  histograms.  Each  of  these  figures  has  three 
parts: 

(a)  the  time  series,  X(t); 

(b)  the  spectral  density,  S(f); 

(c)  the  normalized  histogram,  h(x). 

For  example.  Figure  8a  shows  the  time  series  of  a  section  of  a 
data  file  having  I.D.  number  C506T4.  The  sampling  rate  is  one  per 
second.  Vertical  scale  calibration  is  not  very  relevant  to  the 
present  study;  so  the  vertical  axis  is  omitted.  Figure  8b  shows  the 
spectral  density  of  the  time  series,  which  is  arrived  at  by  computing 
the  FFT  of  the  record  and  plotting  its  magnitude  out  to  0.5  Hz.  The 
spectral  density  of  record  C506T4  is  unimodal  because,  if  one  ignores 
the  fine  structure  of  narrow  peaks  and  valleys,  the  graph  shows  a 
single  "bump’1  spreading  almost  symmetrically  about  0.16  Hz.  Strictly 
speaking,  the  magnitude  of  the  FFT  should  be  squared  to  obtain  the 
spectral  density  as  that  term  is  defined  in  most  texts  on  random 
processes.  Squaring  would  obviously  tend  to  highlight  both  the 
midband  and  the  fine  structure.  Figure  8c  (solid  line)  is  the 
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FIGURE  38.  SPECTRAL  DENSITY  S(f)  FOR  DATA  FILE  P577R6 


FIGURE  3C.  NORMALIZED  HISTOGRAM  h(x)  FOR  DATA  FILE  P577R6 
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FIGURE  5B.  SPECTRAL  DENSITY  S(f)  FOR  DATA  FILE  P340R3 


FIGURE  5C.  NORMALIZED  HISTOGRAM  h(x)  FOR  DATA  FILE  P340R3 
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FIGURE  8C.  NORMALIZED  HISTOGRAM  h(x)  FOR  DATA  FILE  C506T4 
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FIGURE  9B.  SPECTRAL  DENSITY  S(f)  FOR  DATA  FILE  C231T4 
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FIGURE  12C.  NORMALIZED  HISTOGRAM  h(x)  FOR  DATA  FILE  C514T4 
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NORMALIZED  HISTOGRAM  h(x)  FOR  DATA  FILE  A508T2 
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normalized  histogram,  or  sample  p.d.f.,  of  the  same  record.  The 
limits  on  the  horizontal  axis  are  determined  by  computing  the  r.m.s. 
value  of  X(t)  and  multiplying  it  by  +5  (for  the  right  hand  limit)  and 
-5  (for  the  left) .  The  dynamic  range  of  the  abcissa  is  then  ten  times 
the  r.m.s.  fluctuation.  This  range  is  subdivided  into  100  "bins"  of 
equal  width,  with  50  bins  on  each  side  of  the  vertical  axis  (x  -  0). 
The  two  or  three  thousand  sample  values  are  sorted  into  the  bins  to 
obtain  a  histogram.  The  dotted  curve  shows  the  Gaussian  p.d.f.  with 
standard  deviation  equal  to  the  r.m.s.  fluctuation. 

Table  1  summarizes  the  authors'  observation  on  the  20  trios  of 
charts  contained  in  Figures  3  through  22.  The  first  two  columns  in 
the  table  establish  the  correspondences  between  the  figure  numbers  and 
the  file  I.D.  numbers.  The  third  column  indicates  the  modality  of  the 
spectral  density.  A  "1"  in  the  third  column  means  there  is  only  one 
distinct  center  frequency.  A  "2"  in  the  third  column  says  the 
spectral  density  is  Bimodal,  as  illustrated  most  clearly  in  Figure  9b, 
where  the  two  modes  appear  to  be  at  around  0.07  Hz  and  0.17  Hz.  If  a 
"1*"  appears  in  the  third  column,  it  refers  to  the  fact  that  the 
spectral  density  is  unimodal  and  asymmetric,  declining  more  sharply  to 
the  left  of  the  peak,  more  gradually  to  the  right.  The  Bretschneider 
spectrum  and  most  of  the  other  theoretical  ocean  wave  spectra  predict 
this  kind  of  unimodal  asymmetry.  The  fourth  column  in  the  table 
comments  on  the  proximity  of  the  normalized  histogram  to  the  Gaussian 
p.d.f.,  with  "1"  indicating  a  reasonable  fit,  and  "0"  indicating  a 
"poor  fit."  These  entries  are  somewhat  subjective  and  the  reader  may 
reach  different  conclusions. 

Regarding  the  origins  of  the  data,  those  files  whose  I.D.  numbers 
begin  with  the  prefix  "P"  were  acquired  in  a  sheltered,  shallow, 
coastal  location.  The  others,  with  prefixes  "C"  and  "A", 
acquired  in  the  open  ocean. 


were 
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TABLE  1.  OCEAN  WAVE  DATA 


Figure  No. 

File  I.D. 

Spectral 

Modality 

Normalit 

3 

P577R6 

1 

1 

4 

P345R4 

2 

0 

5 

P340R3 

1 

1 

6 

P584R6 

1* 

1 

7 

P584R4 

1 

1 

8 

C506T4 

1 

1 

9 

C231T4 

2 

(0) 

10 

C540T4 

(2) 

(1) 

11 

C516T4 

1* 

1 

12 

C514T4 

1 

0 

13 

C237T4 

(1*) 

(0) 

14 

C211T4 

(2) 

1 

15 

C183T4 

1 

(1) 

16 

C558T4 

1* 

1 

17 

C518T4 

1 

0 

18 

A633T4 

1 

(1) 

19 

A455T4 

1* 

1 

20 

A515T2 

1 

0 

21 

A508T2 

1* 

0 

22 

A380T2 

1 

0 
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NON-GUASSIAN  FLUCTUATIONS 


A  variety  of  techniques  have  been  developed  by  statisticians  for 
measuring  the  deviation  of  a  sample  distribution  from  the  Gaussian 
"standard."  Special  graph  papers  are  commonly  employed.  The  analyst 
may  compile  the  distribution  of  the  squared  sample  values  to  perform  a 
chi-squared  goodness-of-f it  test,  since  the  squared  Gaussian  random 
variable  has  the  chi-squared  distribution.  Another  technique,  which 
derives  from  communications- theoretic  considerations,  is  to  plot  the 
function 


Q(x)  =  -  df  lo9  P(x)  (14 

where  p(x)  is  the  best  estimate  of  the  p.d.f.  of  x,  given  the  sampled 
values  of  X.  For  if  p(x)  is  Gaussian,  one  has 


d  T  1  2,„_2l 

'  dx  log  [  V2TT  e_X  /2°  J  =  x/°2  r 

_  9 

a  straight  line  through  the  origin  with  slope  of  a  .  Thus  the 
proximity  of  Q(x)  to  a  straight  line  with  slope  equal  to  the 
reciprocal  mean  square  fluctuation,  reflects  the  closeness  of  the 

O 

distribution  to  the  Gaussian. 


The  determination  of  the  "best  estimate,"  p(x),  from  the 

normalized  histogram,  h(x),  is  the  subject  of  p.d.f.  estimation 

theory.  Perhaps  the  simplest  procedure,  which  applies  to  cases  in 

which  the  number  of  sample  values  far  exceeds  the  number  of  "bins" 

used  to  compile  the  histogram,  is  distr ibution- independent  linear 
9 

Let  the  bins  be  centered  on  the  regular  sequence  of 


smoothing . 
abcissas  X 


1'  X2 


•  •  r  X  j  , 


.,  X2jf  spanning  the  full  dynamic 
range,  with  Xj  +  xj+i  =  0*  Let  Nj  be  the  number  of  sample 
values  of  X ( t )  which  fall  into  the  j-th  bin.  Then  if 
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2J 

N  -  Z  V 

j=i 

the  normalized  histogram  is  the  sequence  of  numbers 

H j  =  N j/N , 

with  Hj  corresponding  to  h(Xj).  To  symmetrize  the  sample  p.d.f., 
let 

Pj  *  °'5  <Hj  +  ^J+l-j1' 


so  that 


Pj  "  P2J+l-j* 

Now  "smooth"  the  symmetric  sample  p.d.f.  using 

+  °-2Pj-l  +  °-4Pj  +  °'2Pj+l  +  0 ' lp j+2 
or  any  similar  rule  that  preserves  the  identity 


2J 


E 

j-i 


=  l 


This  procedure  was  applied  to  six  of  the  normalized  histograms 
shown  previously;  and  the  results  are  presented  in  Figures  23  through 
28.  These  are  the  smoothed  sample  p.d.f.’s  for  data  files  A515T2, 
A508T2 ,  A380T2 ,  A455T4,  C514T4,  and  P345R4,  respectively. 


The  final  step  is  to  graph  the  six  corresponding  functions  Qj 
which  represent  Q(x),  defined  in  equation  (14).  Because 


d  (logf ) /dx  =  (1/f)  df/dx, 


SMOOTH  SAMPLE  PDF  FOR  DATA  FILE  A380T2 


FIGURE  27.  SMOOTH  SAMPLE  PDF  FOR  DATA  FILE  C514T4 


FIGURE  28.  SMOOTH  SAMPLE  PDF  FOR  DATA  FILE  P345R4 
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one  has 


0^  =  (2/D)  -  PjJ/tPj.i  +  Pj) 

in  the  finite  difference  formalism,  where  D  is  the  width  of  each  bin 
in  the  original  histogram.  These  results  are  shown  in  Figures  29 
through  34,  with  Figure  29  corresponding  to  Figure  23,  and  so  on.  In 
each  case,  there  is  considerable  "choppiness"  in  Q(x)  for  x  near  the 
ends  of  the  horizontal  scale.  The  most  nearly  Gaussian  of  the 
smoothed  sample  p.d.f's  appears  to  be  that  corresponding  to  data  file 
A455T4  (Figure  26).  It  is  not  surprising  that  Q(x)  in  Figure  32, 
corresponding  to  the  same  file,  follows  the  straight  line  rather 
closely  over  the  central  portion  of  the  range.  Figure  23, 
corresponding  to  data  file  A515T2,  appears  more  compressed  in  the 
center,  and  more  persistent  in  the  tails,  than  the  Gaussian  density 
(indicated  by  the  discrete  x-marks) .  Accordingly,  the  plot  of  Q(x)  in 
Figure  29  shows  an  S-curve  in  the  central  range  and  a  clockwise 
curvature  through  the  straight  line  outside  the  range.  It  can  be 
shown  that  the  clockwise  bending  of  Q(x)  symptomizes  "impulsiveness" 
in  the  underlying  random  process,  while  counter-clockwise  bending 
means  the  fluctuations  are  more  confined  and  that  large  fluctuations 
are  less  common  than  for  a  Gaussian  process. 


CONCLUSION 


The  data  presented  in  this  publication,  which  were  not 
pre-selected  or  screened,  provide  limited  support  of  ocean  wave 
theories  that  predict  a  Gaussian  distribution  of  fluctuations.  On 
the  other  hand,  those  data  which  were  subjected  to  a  more  detailed 
examination  in  the  preceding  section  show  a  higher  probability  of 
large  fluctuations  (several  times  the  r.m.s.)  than  predicted  by  a 
Gaussian  model.  The  data  files  in  this  latter  set  were  chosen  to 
highlight  the  disparity. 
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Large  fluctuations  are  most  interesting  to  the  engineer  or  analyst 
who  is  concerned  with  the  design  of  an  oceanic  structure  or  system 
possessing  a  sensitivity  to  waves  that  is  an  increasing,  nonlinear 
function  of  the  amplitude.  In  such  cases  as  these,  one  would  rather 
have  a  p.d.f.  that  yields  a  good  approximation  to  the  extreme  value 
distribution,  even  if  the  matchup  is  poor  at  low  amplitudes.  While 
the  data  presented  here  could  be  used  to  estimate  the  extreme  value 
distribution,  confidence  in  the  result  would  be  rather  low  owing  to 
the  insufficiency  of  the  volume  of  data.  Indeed,  19  records  of  15 
minutes  each  add  up  to  4.75  hours  of  observation;  whereas  the  engineer 
may  be  concerned  with  systems  that  experience  fluctuations  for  weeks 
or  years  consecutively. 
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